Control of Material Damping in High-Q Membrane Microresonators 
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We study the mechanical quality factors of bilayer aluminum/silicon-nitride membranes. By coating 
ultrahigh-Q Si3N4 membranes with a more lossy metal, we can precisely measure the effect of material loss 
on Q's of tensioned resonator modes over a large range of frequencies. We develop a theoretical model that 
interprets our results and predicts the damping can be reduced significantly by patterning the metal film. Using 
such patterning, we fabricate Al-SisN4 membranes with ultrahigh Q at room temperature. Our work elucidates 
the role of material loss in the Q of membrane resonators and informs the design of hybrid mechanical 
oscillators for optical-electrical-mechanical quantum interfaces. 
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Silicon nitride membranes have recently emerged as 
promising resonators for applications ranging from precision 
sensing to realization of a mesoscopic quantum harmonic os- 
cillator OJ-fJl]. Because of their large tensile stress, SiaN 4 
membranes can have MHz resonant frequencies with sub-Hz 
damping rates. The resulting room-temperature Q-frequency 
products of above 10 13 Hz approach the performance of 
quartz oscillators J3|. This ultrahigh Q combined with a two- 
dimensional geometry is an ideal platform for control and de- 
tection of motion in a high-finesse Fabry-Perot cavity, and 
cooling SiaN4 membranes in such cavities to their quantum- 
mechanical ground state is a near-term prospect. However, 
the mechanisms that limit the realized Q-factors of these ten- 
sioned resonators are just beginning to be explored §^f^. 

To date, studies have focused on pure dielectric Si3N 4 
membranes, but a variety of proposed cavity mechanics exper- 
iments would be enabled by the addition of a metallic layer to 
SiaN4 while maintaining high Q &I2]. Foremost, a metallic 
membrane section could form a capacitor plate that couples 
to a microwave LC resonator; in fact, pure metallic drums 
have recently been ground-state cooled using a combination 
of cryogenic and microwave cavity cooling |9|]. With a hy- 
brid dielectric/metallic membrane, one could couple mechan- 
ical motion simultaneously to optical light and a microwave 
electrical circuit in the quantum regime fioll . Such a device 
could solve the difficult, yet crucial, problem of transferring 
quantum states between microwave and optical photons. It 
could also enable enhanced detection of excitations in a room- 



temperature electrical circuit via photodetection 111 111 . Further, 
a magnetic metallic film could be used for magnetic coupling 
of spins to membrane motion 11211 . However, the success of 
these applications will hinge on creating hybrid membranes 
with a sufficiently high quality factor at relevant temperatures. 

In our work, we add metallic thin films to Si3N 4 mem- 
branes and explore the quality factor of many spatial modes 
of the membranes (Fig.Q]). First, we identify two distinct loss 
mechanisms in our experiments: (1) loss of energy from the 
mechanical mode into the substrate, i.e. radiation loss, and (2) 
material damping due to the lossy metallic film. Then, we iso- 
late the material loss-limited Q and develop an anelastic the- 
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FIG. 1: (color online). Geometry of membrane modes, (a) Image 
of a 1 mm membrane in its silicon frame, (b) Illustrations of the 
(m, n) = (2,2) and (15, 15) modes, (c) Image of a patterned Al 
film on top of SiaN4. The central square is Al and the rim is Si3N4 
suspended on a Si frame, (d) Schematic diagram of the Al/Si3N4 
bilayer membrane (50 nm Al and 50 nm SisN4). 



ory that explains the observed dependence of Q on frequency 
for a general clamped, lossy membrane. Our work clarifies 
the role of material loss in highly-stressed two-dimensional 
resonators, and has significant predictive power. Finally, we 
calculate and demonstrate that by removing the metal in a very 
small region near the clamp, we can create metallic SisN4 
membranes with impressive quality factors of over 5 x 10 6 at 
1 MHz at room temperature. 

We use 50 nm-thick stoichiometric LPCVD nitride mem- 
branes that are supported by a 200 |im-thick silicon frame 
(from Norcada Inc.). The membranes are in a square geom- 
etry of side length I = 0.5 mm or 1 mm with tensile stress 
a ~ 0.9 GPa and mass density p ~ 2.7 g/cm 3 . The mem- 
brane mode shapes are given by approximately sinusoidal 
functions like those shown in Fig. [Tib) with resonant frequen- 
cies f mn ~ y/a(m 2 + n 2 )/4pl 2 , where to, n are the integer 
mode indices representing the number of antinodes. The sili- 
con frame is glued at three corners to a metal form on a piezo- 
electric actuator. To probe the mechanical displacement, we 
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have higher Q than asymmetric modes, as predicted by a radiation f (MHz) 

loss model @]. 



position the membrane at the end of one arm of a Michelson 
interferometer. We characterize the mechanical quality factor 
by monitoring the ringdown of the mechanical excitation as a 
function of time in vacuum of less than 10 -6 torr. 

In a first experiment, we measured the Q's of pure SiaN4 
membranes. As shown in Fig. [2] we have the ability to mea- 
sure the quality factors of many modes (up to 150) with differ- 
ent symmetries and to confidently assign a mode (m, n) to all 
measured points. When the data are plotted versus resonant 
frequency (green circles in Fig. [3), the Q is non-monotonic. 
However, when the data are plotted as a function of mode in- 
dex in each dimension (Fig. |2), we see that the asymmetric 
modes (indices n and m dissimilar) have strikingly smaller 
quality factors than do the more symmetric modes (n and m 
nearly equal). This observation is consistent with expected 
trends for radiation loss of elastic waves through the mem- 
brane clamp. As recently calculated and measured in Ref. JH, 
the degree of destructive interference of elastic waves in the 
substrate is responsible for the symmetry dependence. While 
we see consistently low Q's with highly asymmetric modes, 
we see some variability due to the membrane mounting struc- 
ture especially among the lowest order modes @], as expected 
for a radiation loss mechanism jjj. However, for the high- 
order symmetric modes that asymptote to Q's over a million, 
it becomes unclear whether radiation or material loss is the 
dominant effect. 

In our next experiments, we deposit 50 nm of Al using e- 
beam evaporation on top of the pure SisN4 membrane mea- 
sured in Fig.|2jb). The membrane remains under large tensile 
stress, but adding the additional film does decrease the effec- 



FIG. 3: Extracting the material loss-limited Q. (a) Measured quality 
factors of a square SisN4 membrane before (green circles) and after 
(blue squares) adding a 50 nm film of Al. The modes limited (not 
limited) by radiation loss are marked by open (closed) circles. The 
data are plotted as a function of frequency measured after adding the 
Al. (b) Mechanical linewidth of the modes limited by the material 
loss of Al for 0.5 x 0.5 mm (red diamonds) and lxl mm (blue 
squares) membranes. To compare to theory, we calculate the damp- 
ing rate 7 mn /2-7r for each mode, and the points are connected with 
the displayed lines. 



tive stress to <7 e a = 0.35 GPa. With the addition of the metal, 
we see a drop in Q to a maximum of ~ 2 x 10 5 as shown in 
Fig- Ha). Again, as a function of frequency, the Q's are non- 
monotonic, but by drawing on our knowledge of the Q's of the 
pure SiaN4 membrane [Fig. 2(b)], we can clearly distinguish 
radiation loss and material loss. The open squares in Fig. [5J a ) 
represent the asymmetric modes found to be radiation-loss 
limited for the pure S13N4 membrane. If we remove these 
points from the Al-SiaN4 membrane measurements, we ar- 
rive at a clean set of points (closed squares) representing the 
material loss-limited Q as a function of frequency. Damping 
rates 7 = 27r/ / Q for two datasets obtained using this method 
are shown in Fig.[3b). 

We have developed a theoretical framework to describe the 
frequency dependence of the material loss-limited quality fac- 
tors of our two-dimensional structures. We model the mem- 
brane as an anelastic plate that dissipates mechanical energy 
under cyclic loading fl 13M . Under oscillation, the material's 
strains and stresses are not in phase, and the energy supplied 
by the out-of-phase stresses is converted irreversibly to heat. 
This picture has been successfully developed to understand 



3 



damping in one dimensional S13N4 strings lfl4l[l5ll . 

In our case, we start by applying standard plate theory 
with an in-plane force lfl6ll . i.e., under tensile stress, to deter- 
mine the normal modes. The modes must satisfy the bound- 
ary conditions of the clamped plate W = (d/dx)W — or 
W~ (d/dy)W=0foi all four edges. We express the 2D mode 
function W mn (x, y) as a product of stressed-beam functions 
u m (x)u n (y)- We have verified the accuracy of this descrip- 
tion via perturbation theory Jl7ll . We use a closed-form ex- 
pression for the function u„(x) that is a sinusoid with an ex- 
ponential correction near the edge for the clamped boundary 
condition. 

For each mode, we can calculate the loss due to anelas- 
ticity. The oscillation of the plate induces oscillating strains 
£xxS zut , e yy e' lu ' t , and e xy e lut , and the accompanying stresses 
are given by the usual constitutive equation of classical 
plate 1 1 811 with the complex Young's modulus E = E\ + iE%, 
where E 2 is called the loss modulus. During one cycle, the 
full expression for the energy lost is 



MJ= 



2TrE 2 (x,y) ( {e^+eyyf 
l + v I 2(1 - v) 



4 1 
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where v is the Poisson's ratio 111 711 . Note, the strain term 
e xx = — z(d 2 W/dx 2 ) is proportional to the curvature of the 
mode function. To calculate the quality factor, we also need an 
expression for the total stored energy. It can be obtained from 
the maximum kinetic energy U = 2pn 2 f 2 J W(x, y) dV . 
The quality factor for a particular mode W mn is then given 
by Qmn = 27r[/ mn /A[/ mn . 

We start by using our theory to calculate the damping of 
fully-metallized membranes. We apply a least-squares fit to 
the two datasets (two different-sized membranes) in Fig. [3b) 
assuming a single frequency-independent loss modulus. This 
reveals an effective bilayer E2 = 0.55 GPa. The correspond- 
ing Al loss modulus is consistent with typical values for thin- 
film polycrystalline Al at room temperature, as measured, for 
example, via depositing Al on a low-loss Si cantilever 11911 . 
The presumed microscopic origin of the loss is related to crys- 
tallographic defects such as grain boundary sliding 1 2^, 21 ] or 
kinks on dislocations 12211 . Despite this underlying complex- 
ity, our model assumes very little about the microscopic origin 
of the loss. Namely, we assume that the defects are uniformly 
distributed within the deposited metal in the x and y direc- 
tions. We also assume the temperature stays sufficiently con- 
stant in our measurements so as not to affect the loss modu- 
lus. We have verified that the heating due to our measurement 
laser of power 150 |J.W is not a significant effect by measur- 
ing constant quality factors as the power is varied from 10 to 
900 |aW. 

With continued analysis of the theory we can not only 
model, but understand the Q dependencies seen in Fig. [3] and 
put our observations in the context of other studies in ID and 
2D EH [3G3]. We would like to understand: (1) The fre- 
quency dependencies, i.e., why an extremely corrugated mode 
has only a slightly lower Q than the fundamental mode in our 



measurements (2) The geometry dependence, i.e., how damp- 
ing should scale with resonator size. First, we address the 
frequency dependence. As noted above in the discussion of 
Eq. (1), the loss is given by an integral of terms proportional to 
the mode curvatures. We identify two contributions to the cur- 
vature, namely that induced at the clamped edge and that near 
the antinodes in the interior of the membrane. If the curvature 
at the edge dominates we expect a fiat Q as a function of fre- 
quency, or if the antinode contribution dominates we expect a 
decreasing Q as the frequency (and correspondingly the num- 
ber of antinodes) increases. We quantify these statements by 
deriving a simplified expression for Q as a function of mode 
indices m and n for the limit of (1) an isotropic membrane, 
i.e. constant E2 in x and y and (2) high-stress quantified by 
small Am and An where A = y E'h 2 /3al 2 is a dimensionless 
stress parameter. Here E' = Ei/(l — v 2 ) and h, /, and a are 
the height, length, and stress of the membrane respectively. In 
these limits, Eq. (1) becomes an integral over squared sinu- 
soidal terms (antinode contribution) and an exponential term 
(edge contribution) to give a total Q of II 1711 



Q, 



XE, 



1 + A^- 
x 4 

edge antinode 



(2) 



The term A(m 2 + n 2 )ir 2 /4 determines whether there will be 
a frequency-dependent Q. For our experiments, and similar 
experiments with large membranes JiHg], A ~ 10~ 4 — 10~ 3 , 
and hence we expect a relatively fiat Q. However, if the edge 
length is decreased, A increases and the antinode contribution 
can become large. Hence, a frequency dependence appears for 
experiments such as those in Ref. Ill4tl where shorter strings 
(< 35 |^m) are used. 

Further, the prefactor 1/ A in Eq. (2) determines the geome- 
try and stress dependence for the Q of the fundamental mode. 
Physically, A can be written as the ratio of bending energy to 
elongation energy lfl7ll . and as discussed in Ref. lll4l . excit- 
ing energy in the form of elongation energy rather than bend- 
ing energy leads to higher Q. More concretely, based upon 
Eq. (2), we predict that if the membrane side length I is dou- 
bled, Q of the fundamental mode will double for the same loss 
modulus, and this is exactly what is observed in Fig.[3jb); the 
analogous effect in ID was observed in Ref. lfl4T . We also 
see that as the stress is varied, Q scales with ^/a, and hence 
the linewidth 7 = 2nf / Q remains constant. While Eq. (2) 
only holds in the stressed limit, a calculation in the zero stress 
(flexural) limit reveals the linewidth increases by a only a few 
factors from the highly-stressed case. 

Our analysis above indicates that by making the loss mod- 
ulus near the membrane edge small, we can reduce the loss 
significantly. Using our ability to control the addition of ma- 
terial loss with the Al film, we can directly test this prediction. 
The inset to Fig. [He) shows a 1 x 1 mm membrane where we 
deposited Al nearly everywhere except in a small ~50 \im re- 
gion near the edge. The quality factors of this membrane were 
measured to be dramatically higher (blue circles) than a con- 
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FIG. 4: (color online). Ultrahigh-Q metal-covered membranes. 
Measured quality factors of two Si3N4 membranes with Al every- 
where but near the edge of the membrane (blue circles) [see Fig.fTJc)] 
and a full film of Al as a control experiment (red squares). Asymmet- 
ric modes (with one mode index less than or equal to two) are marked 
by open circles or open squares. Calculated quality factors for each 
geometry are shown by the two lines; a single loss modulus is used 
for both. (Inset) The square of the curvature of a stressed mode as 
a function of distance along one coordinate of the membrane. This 
function decays exponentially near the membrane edge. 



trol experiment (red squares) in which an identical layer of Al 
was deposited everywhere on a separate membrane (Fig. 0J. 
For these data we show the Q's measured for all modes, but 
identify the lower-Q asymmetric modes by open circles or 
squares. For both datasets in Fig. [4] we anneal the membranes 
at 340°C after depositing the Al film resulting in an effective 
stress of a c s = 0.6 GPa. While annealing was not necessary 
for studying the fully-metallized membranes of Fig. [3] the un- 
equal stress of the Al film in the partially-metallized mem- 
brane makes the higher-order modes difficult to identify. The 
annealing mitigates this problem, but we still cannot identify 
modes past 3 MHz. Hence past this point we analyze the Q 
envelope by measuring the highest Q mode in every 50 kHz 
window. 

We can again apply our theory to quantitatively predict the 
Q for this new geometry. Since Ei now becomes a function 
of position on the membrane, we return to using the full ex- 
pression of Eq. (1). We assign a finite loss modulus for the 
Al region and zero loss for the SiaN4 rim. We approximate 
the mode functions as those expected for a uniformly stressed 
membrane. This reveals the solid line in Fig.@] which we find 
scales as 1/n 2 for diagonal modes, as expected from antinode 
contributions. The dashed line shows the corresponding result 
for the fully metallized membrane using the same loss modu- 
lus. (The loss modulus found here is E<z = 0.3 GPa, which is 
a smaller value than for Fig. [3] due to the annealing.) We find 
the theory successfully traces out the envelope of the mea- 
sured Q's. The lowest order modes of the partially-metallized 
membrane reach as high as Q = 6.5 x 10 6 ; this falls short 
of the predicted Q just from Al material loss (blue line) likely 



because these modes are now again limited by radiation loss. 
In comparison, recent measurements of metal microstrings at 
room temperature revealed Q's of 10 3 — 10 5 J23,24], and even 
at cryogenic temperatures, where the metal's material loss is 
significantly reduced, observed Q values for tensioned mi- 
croresonators are typically 10 5 — 10 6 125142711 . 

It is elucidating to understand what would happen to the Q 
trends we observe for the partially-metallized membrane upon 
varying the membrane stress. This requires analysis of the 
spatial dependence of the curvature in the membrane plane. 
In the inset to Fig. |4] we see the high-curvature area only oc- 
cupies a very small ~ 1 region near the edge of a high- 
tension membrane (here we use our lowest a c s = 0.35 GPa); 
specifically, the decay length is AZ/4 §Y~h . As the stress is 
reduced (and hence A is larger), the curvature becomes more 
uniformly distributed over the membrane plane. Hence we 
would not expect a dramatic difference in Q for a purely flex- 
ural mode when avoiding lossy material only at the edge. 

The localized curvature of the tensioned membrane that we 
observe provides insight into a variety of membrane appli- 
cations. Note for higher-order two-dimensional modes the 
curvature varies along the edge of the membrane, i.e. there 
are low-curvature regions near the nodes at the membrane 
edge (Hi. Thus, to create an electrical link between a cen- 
tral metallized patch and external circuits, and maintain high- 
Q performance, one could tailor metal connections to match 
up with the low-curvature regions near nodes at the mem- 
brane edge lfl7ll . Further, membrane patterning via holes is a 
promising technique to increase reflectivity of membranes for 
optomechanics experiments, but like metal deposition, also 
has potential to introduce defects. A full understanding of the 
curvature of the two-dimensional membrane is important for 
understanding the change in Q, or a lack of a decrease in Q, 
in recent patterning experiments |28, 29ll . 
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SUPPLEMENTARY INFORMATION 
NORMAL MODES OF A HIGH-TENSION PLATE 



For an isotropic square plate with uniform boundary tension and under sinusoidal oscillation in the transverse direction, the 
dimensionless differential equation for a mode W(x, y) is [1] 



_D_(d^W_ d^W_ 



d 2 W d 2 W 



2D d 4 W 

perturbation term 



n 2 w 



(si) 



where £ = x/l, rj = y/l, D is the fiexural rigidity defined by D = Eh 3 /[12(1 — v 2 )\, E is the Young's modulus, h is the 
plate thickness, v is the Poisson's ratio, a is the inplane stress, I is the length of the side of the plate, = uj^J pi 2 /a is the 
dimensionless angular frequency, and p is the mass per unit volume. 

The modes for our device must satisfy the boundary conditions of a plate clamped on all four sides (CCCC), i.e., 



all edges 



dn 



all edges 



= o 



(S2) 



Here n is the axis perpendicular to the edge. 

To solve Eq. dSU and Eq. dS2l ), note that the coefficient of the cross term is small for our devices: 2D j (ahl 2 ) ~ 10~ 6 <C 1. 
Therefore, we can apply standard perturbation theory with e = 2D /(ahl 2 ) as the "small parameter" and d A /d^ 2 drj 2 as the 
perturbation operator. Then the unperturbed equation can be solved by separation of variables. 

\j*W{l-,ri)=X(Z)Y(ri), n 2 



il 2 + Q 2 , and E' = E/(l - v 2 ). We have 



E'h 2 d 4 X 
\2al 2 ~d^ 
E'h 2 d 4 Y 
YLol 2 Ihf ~ ~£hf 



d 2 X _ 2 
d 2 Y 



n 2 Y, 



x 

Y 



77=0,1 



dX 
BY 



€=o,i = 
17=0,1 = 



(S3) 



These are the dimensionless equations of doubly clamped beams with Young's modulus £"and pre-stress a . The eigensolutions 
Um(0 (or u n (j])) and the corresponding frequencies TO (or VL n ) in the limit of small e are given by 



MO. < e < i 

_(-ir+Vn(l-a 5<£<1 



"m(0 



MO = v2a{ sin[/3 m £] H (exp[-a m £] - cos[/3 m £]) 

Q m ~ tyvk 



(S4) 

(S5) 
(S6) 



In Eq. dS5l l a <C lis the vibration amplitude, and a m and fi m are defined as 



1/2 



1/2 



1+ sfl + 2eft= 



1/2 



-1+ y/l + 2eSl* 



1/2 



(S7) 



These functions have a sinusoidal shape with an exponential correction near the edge for the clamped boundary condition. Beam 
functions when e is intermediate in size are given in Ref. [2]. 

The products {u m (^)u n (r})} are then the solutions of the unperturbed equation. The normal modes of the plate, however, 
must be symmetrized because of the square geometry of the plate 



75' 



" in (0 u n(?7) + u n (£)u m (vi) ) , m>n 

(u m (£,) u n (V) - u n (QUm (v)h m < n 

u m (£,)u m (r)), m = n 

n 2 



(S8) 
(S9) 
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The orthogonality condition is 



o Jo 



(S10) 



Perturbation Correction to First Order 



Now we will calculate the first order correction due to the perturbation term. Note that for m ^ n, Qmn = £lnm, i- e -> the 
modes Wm«(£, T}) and Wnm(£, v) aie degenerate. Applying the degenerate perturbation theory, it turns out that the perturbation 
operator is diagonalized in the basis set {Wmn(£, ">])}■ Thus, we can easily calculate the eigenfrequencies to first order m „ 
and the mode shapes with the first order correction W mn (£, rj). Assuming j3 p /a v <C 1 for p <S {m, n, m', n'}, or equivalently 

e7r 2 max{m, m', n, n'} 2 /2 -C 1, we have 



O( ) 2 



(0) 2 



mi i 
1 pi 



< 1 



< 



2TT 2 e 3 / 2 mm'nn' 



<1, 



Therefore, we can use rj)} as the mode functions for our anelasticity calculation. 



(Sll) 



Mode Functions of Bilayer Plate 

For the fully or near fully metallized plates, we calculate mode functions by using a single-layer plate with thickness 100 nm, 
an effective stress, Young's modulus, mass density, and Poisson's ratio. The effective Young's modulus E = E\ = 135 GPa is 
the average of the Young's modulus of Si,3N4 and Al weighted by their thickness. The mass density p of Al and S13N4 are very 
similar, and hence we use an effective mass density of p ~ 2.7 g/cm 3 . We use an effective Poisson's ratio of v = 0.31. 



ANELASTIC LOSS OF A HIGH-TENSION PLATE 



We place the xy plane at middle plane of the plate and let the plate oscillate with a small amplitude at an angular frequency lu 



in the z-direction. The oscillation of plate We 1 
oscillating strains of 



(t) 



induces a variation of local bending and overall elongation [3], and thus gives 
/ -z- 



- !)!) 



Sxy(t) 



\ 



■ dy 

-2z^ 



dxdy 



We 



1^; 



0,1 



(S12) 







local bending-o(W) elongation~e>(w 2 ) 
/ -i 



\ 



, d 2 W 
'~dx r 
„ d 2 W 

W 

2 z d w 

dxdy 




e lut = e Q e lult 



(S13) 



The strain variation is mostly induced by pure bending, as shown in Eq. ( IS 1 2b . The accompanying variation of stresses are given 
by the usual constitutive equation of classical plate [1] with the complex Young's modulus E = E\ + iE 2 

a xx (t) \ E ( 1 v 

*„„(*) = — - v 1 I e(t) = EMe(t) (S14) 
Txylt) ) V V (1 - v)j2 



*(<) = 



Note that the stresses can be separated into in-phase and out-of-phase terms as follows 

»[*(*)] = flt[EMe(t)] = ^[EMe e iwt ] 

= EiMe cos(wi) - E 2 Me sin(wt) 



(S15) 



in-phase 



out-of-phase 



8 



The mechanical work done per unit volume per oscillation is then given by: 

Aw = j^[a T ]^[^\dt 

2tt 

= I [Ei cos(wt)e T M - E 2 sin(o;i)eo T M] sR[e 
Jo 

= TrE 2 e T Me > 



de* 



dt 



-}dt 



This leads to a dissipation. Meanwhile, the bending energy supplied by the in-phase stresses is stored in the system. 



^bending 



U[a(t)} ■ R[^]dt = E 1 s t Me /2 = §-^ 



(S16) 



(S17) 



Eq. (IS 17b is consistent with the strain energy density of a pure bending plate derived in Ref. [4]. We can see that the local 
dissipation Aw is proportional to the density of stored bending energy. In other words, the bending energy supplied by the 
out-of -phase stresses is converted irreversibly to heat. 

To find the total dissipation per cycle, we integrate Eq. dS 16b over entire volume of the plate and insert the strain-displacement 



relations [Eq. (IS12H . We have 



AU 



AwdV -- 

nE 2 (x,y) 
(1-v 2 ) 



TTE 2 eo T Me dV 

{(£ XX + Eyyf - 2(1 - v)(e xx £ vy - e 2 xy / 4)} dxdydz 



2ttE 2 (x, V ) r(e 
[\ + v) I 2(1 - v) 



e 2 i 



z 2 dz 



TrE 2 (x,y) J f cPW_ ^ d 2 W 
dy 2 



dx 2 



-2(1-1/) 



(l-i/*) ^ 

(mean curvature) 2 

The total loss is given by the integration of curvatures over the plate plane. 



dx 2 dy 2 \ dxdy J 



>dxdy 



(SI 8) 
(S19) 

(S20) 



Gaussian curvature 



Stored Energy 

In order to evaluate the quality factor Q = 2ttU/AU, we need to calculate the stored energy U and energy loss AU. We 
calculated AU above, and here we analyze the stored energy, which can be obtained from the maximum kinetic energy or the 
maximum elastic energy. The maximum kinetic energy is 

[Wic = J J W(x, yfdxdy = 2 P ir 2 f 2 J W(x, yfdV (S21) 

The maximum elastic energy for a tensioned resonator is the sum of the displacement-induced elongation energy and the bending 
energy [3]. These two energies are given by 



f/olongation =—JJi j + ) } dxdy 



ah rr fdwy few 

[/bending = J U hcnd ingdV = J ^-^~ dV ( S23 ) 



Effective Loss Modulus 



Using the expressions for U mn and AU mn above, and assuming E 2 (x, y) is constant over the metallized area, we can evaluate 
the quality factor or damping rate for each mode with only one unknown quantity E 2 . In our data analysis we extract an effective 
E 2 for the 100 nm bilayer membrane using a least-squares fit to our data, as discussed in the main text. The loss modulus of Al 
can then be estimated by E 2 (ha,- l . 3 ^ i + ft-Ai)/^Ai- Where /isi 3 N 4 and /jai are the thickness of the SiaN4 and the Al, respectively. 
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INTERPRETATION 
Small Am (An) Limit 



Define a dimensionless parameter A = V2e = \J E' h 2 / ial 2 . In the limit of small Am, which holds for our membranes and 
other high-tension resonators [3, 5-7], Eq. dS7b becomes 



1/2 



1/2 



1+ s/l + 2efl 2 



1/2 



1/2 



-1 + V / l + 2e0| 

The beam equations Eqs. dS4b - dS6b become (back to coordinate .t) 

v m (x), 



1/2 



u m (x) 

t) m (i) = y/2a I sin 



< x < \ 
i 



{~l) m+1 v m (l~x), l<x<l 
mnx 



I 



2 

Xttitt 



exp 



AZ/2 



cos 



TO7T / (7 



(S24) 

(S25) 

(S26) 
(S27) 



Isotropic Plates 

For a plate with uniform thickness and constant E^, we can simplify Eq. dS20t using Green's theorem. The integration of 
Gaussian curvature is zero for a clamped plate (CCCC). 

ffld 2 Wd 2 W fd 2 W\ T 
J J dx 2 dy 2 \dxdy J 

Thus Eq. ( lS20b becomes, 

AU =m^)JJ[^ + iw) dxdy (S28) 

Note d 2 W I dx 2 is the curvature of the mode in the x direction, and d 2 W/dy 2 is the curvature in the y direction. The sum of the 
two curvatures is the mean curvature. 



r d 2 W dW , / d 2 W dW , 
dxdy = f -^r-^-dy - f -^-^--^—dx = 
J oy z ox J oxoy ox 



Curvature Induced near the Clamped Edge and around the Antinodes 

In our anelastic model, the local dissipation of a clamped plate is proportional to square of the mean curvature, as shown in 
Eq. (lS281 l. Hence, it is important to understand its spatial dependence. Since our 2D normal modes are the product of ID modes, 
we start by calculating the curvature of a ID mode u n . In the limit of Art <C 1, we can approximate the square of the curvature 

as 



d 2 u n (x) 
dx 2 



2n 2 TT 2 a ( 4 



I 4 



^exp 



M/4 



n 2 TT 2 sin 2 



/ 



I 

< x < - 
~ ~ 2 



(S29) 



The equation above describes the x dependence of the curvature in the x direction. Note that the coefficient of the first term 
is much larger than that of the second term. This first exponential term is concentrated in a small range XI/ 4 = h^/E' /48a, 
leading to significant loss near the edge. The second term distributes over the entire length with sinusoidal dependence, just like 
the curvature of a simply-supported beam. 

For the 2D case, in addition to the above-mentioned features, the large curvature induced near the edge oscillates sinusoidally 
along the plate edge. To see that, we calculate the mean curvature square of the diagonal mode (n, n) at the edge x = 



10 



d 2 W nn | d 2 w nn 
dx 2 dy 2 



x=0 



16n 2 ir 2 a 2 
X 2 P 



■ sin 



nny 



(S30) 



See Fig.[ST]for an illustration of square of the mean curvature over the entire plate plane. 



Edge Loss and the Frequency Dependence of the Quality Factor 



To understand the role of loss induced near the clamped edge, here we study Q of a mode (m, n) in the limit of small An and 
Am. First, we calculate the maximum kinetic energy. Applying Eq. dS 10b to Eq. dS21l ), we have 



pha 2 l 2 (Lu 2 m + <)/2 - aha 2 {m 2 + h 2 )tt 2 /2 



(S31) 



Second, we calculate the loss per cycle drawing on our earlier curvature analysis. Here we use W mn (x, y) = u m {x)u n {y). We 
have verified that the result is the same as using the symmetrized mode functions in Eq. 

2 



AU„ 



ttE 2 
1-v 2 



-dz 



d 2 w mn , d 2 w mn 



dx 2 



dy 2 



dxdy 



12(1 - v 2 ) 

nE 2 h 3 
12(1 -v 2 ) 



nE 2 h 3 
12(1 - V 2 ) 



d 2 u m (x) d 2 u n (y) 
u n {y) + u m [x) 



2al 



dx 2 

d 2 u m (x) 
dx 2 

d 2 u m {x) 
dx 2 

l ' 2 fd 2 u m {xY 



dy 2 



dxdy 



dx / u n (y) 2 dy+ / u m (x) 2 dx 



u n (y)dy 



d 2 u n (y) 
dy 2 



dy 



u m (x)dx 



d 2 u n {y) 
dy 2 



dx 2 

1/2 d 2 u m (x) 
dx 2 



dx + 2al 



l ' 2 fd 2 u n (y) 
dy 2 



dy 



u m (x)dx 



l ' 2 d 2 u n (y) ] 
— 2 — u n (y)dy 



nE 2 h 3 f 16m 2 TT 2 a 2 f 1 ' 2 



12(l-i/ 2 ) 



PX 2 



cxp 



AZ/4 



dx 



dy 2 

16n 2 ir 2 a 2 f 1 ' 2 
PX 2 



exp 



XI/ 4 



dy 



"edge loss" 



Am it a 
P 



A/2 


'mirx' 


/ sin 2 
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- 1 \ 



dx 



AnWa 2 r 1 ' 2 



P 



rmy 



dy 



2rrm 2 a f l/2 



I 2 



"antinode loss" 



dx 



I 



2rm 2 a f 1 ' 2 



P 



nixy 



dy 



PX 



TiE 2 h 3 
12(1 - I/ 2 

edge 

(m 2 + n 2 )TT 3 E 2 h 3 a 2 
3(1 — i/ 2 )Z 2 A 



"antinode loss r 
4(m 2 + n 2 )7r 2 a2 (m 2 + n 2 ) 2 7T 4 a 2 



P 



1 +A 



antinode 

(m 2 + n 2 )-K 2 



(S32) 



(S33) 



edge antinode 

We see that the ratio of the antinode loss to the edge loss is given by A(m 2 + n 2 )7r 2 /4. Combining Eq. ( 1S3U and dS33K we 
obtain 



Q, 



2wU„ 



AU„ 



3(l-i/ 2 W 2 A/ (m 2 + n 2 )7T 2 

eJ ' 1 + a 



IEi 

XE 2 



1 +A 



(m + n )ir 



(S34) 



edge antinode 
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x (mm) 



FIG. SI: The spatial dependence of the mean curvature square of mode (2, 2) over entire plate plane. A low stress of 0.1 MPa is used and the 
values of curvature are scaled nonlinearly to make the structure more visible. 



Similarly, we can derive the formula for a ID string (with v = 0) 

1 + A—— I (S35) 



\E l (. ,n 2 7r 2N " 



A E- 



For small A as low as 10~ 4 — 10~ 3 in our devices and Ref. [5-7], Eq. ( !S34t and Eq. dS35l ) correspond to a slightly decreasing 
Q as a function of frequency. For A~ 2 x 10~ 2 as in Ref. [3], the edge loss and the antinode loss are comparable when n = 4. 
This gives a stronger dependence on n, and hence frequency. 



Physical Meaning of A 



To investigate the physical meaning of A, we calculate the elongation energy of (m, n) mode in the limit of small An: 



Uc 



loneation 



2 

= ahla 

2ir 2 aha 2 
I 



du m (x) 



dx 

l ' 2 Uv m { x y 

o v dx 

1/2 



Uniy) 2 + U m (x) 



du n (y) 



dx 



dy 

t/2 fdv n (yy 2 

dy 



dxdy 



dy 



I 



dx + n~ 



1/2 



nny 



dy 



= IT aha (TO + n )/2 = t/kmctic = ^elastic 

This means that the maximum elastic energy is dominated by the elongation energy. Thus, 

q ^7rt/ mn 27r[/ c i on g a ti on Ei ^elongation 

AU mn 2-k{EiIE\) ^bonding E 2 C/bending 



Comparing with Eq. dS34b . we find 



t^bendi: 



^elongation 1 

_ A 
1 

~ A ' 



1 + A 



'm 2 + n 2 )-K 2 



for (n, to) = (1,1) 



Thus, A is the ratio of the bending energy to the elongation energy for the fundamental mode. 



(S36) 



(S37) 



(S38) 



REFERENCES 



[1] A. W. Leissa, Vibration of Plates (NASA, Washington, D.C., 1969). 



12 



[2] A. Bokaian, J. Sound and Vib. 142, 481 (1990). 

[3] Q. P. Unterreithmeier, T. Faust, and J. P. Kotthaus, Phys. Rev. Lett. 105, 027205 (2010). 

[4] S. Timoshenko, Vibration Problems in Engineering (D. Van Nostrand Company, Inc., New York, 1937). 

[5] I. Wilson-Rae et al, Phys. Rev. Lett. 106, 047205 (201 1). 

[6] A. Jockel et al, Appl. Phys. Lett. 99, 143109 (201 1). 

[7] S. Schmid, K. D. Jensen, K. H. Nielsen, and A. Boisen, Phys. Rev. B 84, 165307 (201 1). 



